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Abstract 

Using the dual string theory, we study a circular baryonic configuration in a wind of 
strongly coupled A/" = 4 Yang- Mills plasma blowing in the plane of the baryon, before 
and after a quark has dissociated from it. A simple enough model that captures many 
interesting features is when there are four quarks in the baryon. As a step towards 
phenomenology, we compare representative dissociated configurations, and make some 
comments about their energetics and other properties. Related results that we find 
include the observation that the screening length formula LgT ~ (1 — 1>^)^/^ obtained 
previously for other color singlet configurations, is robust for circular baryons as well. 
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1 Introduction and Conclusion 



There is a possibility that screening of heavy quark baryons in a wind of strongly 
coupled plasma might be experimentally accessible at RHIC or LHC0. Unfortunately, 
this is a deeply non-perturbative scenario in standard QCD, and therefore essentially 
out of theoretical control. Interestingly enough, semi- quantitative features of baryon 

^We emphasize that so far heavy quark baryons have not been observed either in elementary 
coUisions or heavy ion coUisions. 
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screening can be computed using a dual string theory through the AdS/CFT corre- 
spondence [Ij. This has resulted in some non-trivial progress in the understanding of 
screening phenomena in strongly coupled plasmas [21 El HI [5] . 

In a wind of plasma, the screening is direction-dependent, and it stands to reason 
that when the baryon dissociates, there will be preferred directions in which this can 
happen first. The purpose of this paper is to explore this possibility by computing 
some basic energetics. This should be taken as a small step towards phenomenology. 

As a warmup, we first compute the (regulated) energy of a circular baryonic 
configuration with Nc = 4: quarkqj (attached in the bulk to a D5-brane baryon vertex 
wrapping the S^), moving in a hot, strongly coupled, Af = 4 plasma. For a generic 
baryon configuration, this can be done using worldsheet string theory in AdS^^ x S^, 
with a black hole in the interior. Our results about these baryons also serve as a test 
of robustness of the known results in the literature. After this, we calculate the energy 
of the configuration after one of the quarks has dissociated from the baryon and is 
well-separated. Both the quark, and what remains of the baryon, will have trailing 
strings reaching down to the black hole horizon. There are many such dissociated 
configurations, and we will try to get some intuition for them by considering some 
special cases. Among these will be two extreme scenarios: one when the remaining 
three quarks are spaced equidistantly in a line along the wind (case I), and the other 
when the quarks (again spaced equidistantly along a hne) are perpendicular to the 
wind (case II). The "dynamics" of the system is sufficiently stringent that for the 
configurations we consider, we will be able to explicitly do the computations without 
getting tangled up in too many details. 

We stress that our aim is not to identify preferred dissociation channels and do any 
detailed phenomenology. This would be a tall order. Our aim is merely to see whether 
something can be said about energetics of transverse configurations vs. longitudinal 
ones. The two specific cases we consider are chosen with this specific purpose in mind, 
and they are not supposed to necessarily be the preferred dissociation products. In 
fact, one might think that the "natural" dissociated configurations are instead the 
ones we consider in Appendix C. But this is so only if one assumes that the remaining 
quarks stay on the circle. The problem here is that "natural" is not very well-defined, 
because the dynamics is not under control and all we are dealing with are static 
configurations. 

^The choice 4 is obtained by optimizing between non-triviality and simplicity. 
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In any event, we find that the energy of the configuration with quarks parallel 
to the wind (case I) is at a higher energy than case II, and therefore it is tempting 
to speculate that for baryons in a wind of plasma, the quarks along the wind will 
dissociate first. We emphasize that the computations we do should not be taken 
as a proof of this claim, even though we believe they are suggestive. In particular, 
it is not clear what configurations constitute extreme cases. It is a complicated 
dynamical question whether a configuration would prefer to change the dimension 
(i.e., change the length L) and/or change the orientation in order to lower its energy. 
But it can still be instructive to have a comparison of identical configurations with 
different orientations: naively, one might think that for a given dimension, a transverse 
configuration is at a higher energy because it has more cross-section to the wind. Our 
results show that this is not the case. Also, there are speculations and comments 
in the literature about the relative energetics of transverse-vs.-longitudinal cases. 
Such statements make most sense only if one assumes that we are comparing two 
configurations of identical dimensions, and this is what we do. 

In general, the features of a dissociated baryon are hkely to depend strongly on the 
specific configurations under consideration, and making generic claims is difficult. To 
emphasize this, we compute the energetics of some other dissociated configurations 
in an appendix. 

A more detailed computation where the number of quarks and the possible dis- 
sociation channels are increased will certainly be useful in shedding more light on 
similar questions. It would be interesting to do a full scan in the the space of allowed 
dissociated configurations, but this is numerically a more demanding problem than 
what we have undertaken in this paper. It should also be pointed out that conclusive 
evidence might require an understanding of the dynamics of dissociation patterns 
which our static calculation is blind to: in principle, it is always possible for instance 
that the configuration can rotate or dilate during or after a dissociation process. 

In the course of the computations in this paper we learn a few things about 
dragging and non-dragging objects and these will be elaborated upon as and when 
they arise. Among these is the observation that the screening length formula found 
previously in the literature for other color singlet configurations, LgT ~ (1 — 
is vahd for circular baryons as well. This is interesting because the details of the 
configurations and the computations are quite different in our case. Another thing 
we find is that the plots of vs. L for non-dragging objects exhibit a cusp while 
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those of dragging configurations exhibit a loop, and we speculate that the area of 
the loop is a measure of the drag of the configuration. This could potentially be a 
measure of the energy loss. 

In this paper, we only consider the simplest possible scenario: a baryon that 
is simulated by a circular configuration of four (external) quarks in the maximally 
symmetric gauge theory, with the added simplification that the plasma is flowing in 
the plane of the baryon. But it is possible to consider generalizations of the results 
here to more generic baryon configurations, less special gauge theories and perhaps 
more generic wind directions. It would be interesting to see how generic the results 
are. Some papers that are relevant (from various angles) to hot strongly coupled 

QCD are[IIl[I21[I3l[ia[I3[ISl[ni[IHl[ISl[2niEIl[221[23|. 

2 Baryon Screening and the AdS Black Hole 

Baryons are made out of fundamental quarks in QCD, but the field content of the 
A/" = 4 supersymmetric Yang-Mills theory consists of gluons, gluinos and scalars (all 
in the adjoint); but not quarks. So in order to model QCD phenomena, we introduce 
baryons that are constructed out of external quarks. If the A/" = 4 SYM theory has 
Nc colors, the baryons will be constructed from such external quarks. 

The dynamics of baryons in the gauge theory is captured in the dual AdS^ x 
string theory through the introduction of the so-called baryon vertex [6J . The claim is 
that baryons in the gauge theory are dual to configurations that involve a D5-brane 
wrapping an in the bulk, with all the heavy external quarks in the boundary 
baryon being linked to it through fundamental strings (all of which are of the same 
orientation). 

The way in which we make predictions for a baryon moving in the plasma is by 
boosting to the rest frame of the baryon and letting the plasma move instead. In the 
dual picture, we look for static baryon configurations in the boosted bulk metric. In 
the course of this paper, we will be exclusively working with the case of A/" = 4 SYM, 
with the plasma at a temperature T. Finite temperature implies that the asymptotics 
of the bulk is still AdS^ x but in the interior we have to change the metric to include 
a black hole whose Hawking temperature is T [?]. Before the boost, this bulk AdS 
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with 



black hole metric takes the form 

ds^ = -f{r)dt' + ^dx' + Tfr + R^dnl (2.1) 

t[r) 

nr) = (1 - ^) . (2.2) 

The asymptotic boundary where the field theory lives is supposed to be at r ^ 00, 
and is spanned by x = {xi, X2, x^} and the bulk-time t. In the above, tq is the black 
hole horizon and the temperature T is fixed by the Hawking relation T = The 
standard CFT to AdS correspondence is usually expressed by taking the defining 
parameters of the gauge theory to be A (= gyM^c) and Nc- Then, the bulk data is 
related to the boundary data through Maldacena's famous relations 

^ = 4„.. VA = :?!. (2.3) 

with the string coupling and the worldsheet tension. The Maldacena conjec- 
ture can be taken as the claim that string propagation with these parameters in the 
AdS background (with a background five-form fiux controlled by R) is just another 
description of the gauge theory. 

We will take the plasma wind to be in the xs-direction after the boost, and the 
velocity and rapidity are related hj v = — tanh?]. The boosted metric takes the form 

ds^ = -Adt^ + 2B dt dx^ + Cdxl + — {dxj + dxj) + j^^dr'^ + R^dilj . (2.4) 

The various quantities are fixed by. 



and 



ri = Tg cosh T], and r2 = Tq sinh rj. (2-6) 



Our description of the baryon-like configurations will involve the baryon vertex 
(the D5-brane), and various strings emanating from the baryon vertex and ending on 
the asymptotic boundary or the black hole horizon. The actions of the string world- 
sheets and the D5-brane can be used to determine static solutions to the equations of 
motion. We will work below in the restricted context of the AdS black hole described 
above, a more general formulation can be found in p]. 
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In this section, we will illustrate the method by explicitly working out the relevant 
physics of an A^^c = 4 baryon moving in the plasma. This example is also our prime 
workhorse. The conflict between setting Nc = 4: and the desire to have a large Nc 
planar approximation where finite string coupling effects are suppressed does not 
seem to be too much, because as we will see, the various results that were found in 
[S] can be reproduced in our case as well. In particular, in [3], large Nc squashed 
baryons were considered, and the curves that we find later in this section are in 
excellent agreement with their worljf]. So we believe that we are not throwing the 
baby out with the bath water by studying this simple system. Adding more quarks is 
in principle straightforward, except that the computational effort is, of course, more. 

The configuration we wish to study is shown in Figure 1. The circular configura- 
tion of quarks lies in the xi-xs plane and the plasma wind is in the xs direction. In 
the figure we have suppressed all coordinates except r,xi and x^. The action for the 
system is given by 

5* = 5*1 + 5*2 + 5*3 + 5*4 + Sd^, (2.7) 

where subscripts denote the various strings and the D5-brane. The actions for the 
various strings are computed a la Nambu-Goto in the black hole background: 

Sng = J dadr^/^g^, with g^p = Gi^^daX^'dpx'' , (2.8) 

where G is the black hole metric and g is the induced worldsheet metric. If we assume 
that the configuration experiences no draj^, we can look for axially symmetric circular 
configurations like the one shown in Figure 1. For the two strings in the Xi direction, 
then, we can take the embedding to be 

r = t, a = r, x'f^ = x^^\a) (2.9) 

where a is either 2 or 4 and denotes the appropriate string (see Figure 1). No r 
dependence arises because we are interested in static configurations. The action for 

■^The details of the configurations we consider here are different from those considered in so 
this resuh also serves as a check of robustness for the screening-length formula. 

This assumption is supported by the computations involving meson configurations [21 [U [H [9l 
[TT1[T2] as well as previously considered baryon configurations [51 [TO]. Our aim in this section is to set 
up the formalism and check whether we can connect with the results obtained in the literature, so 
we will not try to derive it ab-initio. We will assume no-drag for the baryon, and be content that the 
results work out precisely as expected. Of course, when we consider the dissociated configurations, 
we will not make the no-drag assumption. 
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Figure 1: The Nc = 4: baryon configuration in the AdS black hole background. The 
baryon vertex is at r = Tg. 

the strings takes the form 

where T can be thought of as the total time which gets divided out in any relevant 
quantity. Primes denote derivatives with r. In the above expression, re is the position 

of the baryon vertex and it can lie anywhere between r = Tq and r = oo, i.e., between 
the boundary and the horizon. The boundary conditions for the string coordinates are 
fixed by the condition that the baryon vertex has coordinates (r = re,xi = OjX^ = 
0). At the boundary, from a glance at the figure, we see that x\ (oo) = —L and 
a;^'^''(oo) = L, where L is the radius of the circle. For the strings in the Xs-directions, 
similarly, we get 




6 = 1,3. (2.11) 
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The boundary conditions at infinity are given by x^^\oo) 
action for tlie D5-brane vertex can be taken as [5] 

N^RT 



L, X 



(3), 



oo 



S. 



D5 



Stto; 



-L. The 



(2.12) 



2.1 Energetics and Screening Lengths 

Tfie equations of motion for the configuration are obtained by varying the total ac- 
tion with respect to the a;(r)'s and also with respect to the location of the D5-brane. 
The details have been worked out in [3] and the result adapted to our case takes the 
following form. (We suppress the string number superscripts a(= 2,4),6(= 1,3) in 
some of the equations below.) 
Strings 1, 3: 



/ n2 



[Xn 



KIA 



/V2 (rV - R^K. 



2\ ' 



x[ = 0. 



The K3 are integration constants. The K3 have to satisfy the condition, 



Ki'^ + K. 



(3) 



0, 



(2.13) 



(2.14) 



as a force balance condition on the D5-brane. So we effectively need to solve only for 
one of the strings. 
Strings 2, 4: 



/ \2 



R' 



[X. 



r2 / {r^A - R'^Kf) ' 



Again there is the relation 

that has to be satisfied. 
D5-brane: 

2R^/A 



Kf^ + i^f ) = 



Xo = 0. 



(2.15) 
(2.16) 



V/ (i?2 + /r2( 



X 



I \2\ 



2RA 



4 I 4 
r2 + 



(2.17) 



,^ V/ (Ai?2 + /2r2(a;'3)2) 

It should be noted that this equation is evaluated at r = rg. Since the string equations 
written down above depend only on the square of the i^'s, and these squares are 
identical for both strings in each pair, we don't specify the superscript in the x\ and 
the X3 in the D5-brane equation. 
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By introducing new variables, we can write these in the equivalent form 



i=^r,,_i_. (2,18) 

Try, "(y'-p')]! yi-p'-fJ2 > ' 

^ ' ^ =, (2.19) 



ttT 



(y^ - P^)i.y^ - cosh^ T] - a2) 



a/1 — cosh^ t^a/I — — /?2 a/1 — — p"^ cosh^ rj 1+p^ cosher] ('2 20) 
(1^7) ^ ~ 2v/l-p^cosh^' 

We have integrated the string equations of motion while imposing the boundary 
condition that the radius of the circle is L at the asymptotic boundary. In the 
process we have also introduced the notation 

«' = ^, f3' = ^, y = -, P = -. 2.21 

To extract information about the baryon we need to solve these three equations 
simultaneously. This can be done numerically. 

• Pick a value for rj first. 

• Now, for each p (it is easy to see that p must lie in the range (0, 1) from the 
geometry) we can solve for (3 in terms of a as we vary a (which again has to be 
between and 1.). 

• The correct value of a for each p, is the one where the a-integral matches the 
/3-integral. 

• The value of the integral(s) at which this match happens is the value of L. 

• Repeat the above procedure for another value of r]. 

Plot of L as a function of p are presented here, for a few values of t]. The plot for 
L has a peak, and this is what one identifies as the screening length, Lg. The results 
here match perfectly with those in [5]. The dependence of Lg on the rapidity t] can 
also been plotted, and the result again agrees with expectations from other baryonic 
and mesonic configurations considered previously in the literature. In particular, it 
is evident from this plot that LgT ~ (1 — f^)^/^ for large boosts, and this is a check 
of the robustness of this result. 
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0.4 



0.6 



Figure 2: LT vs. p for r/ = (top curve), 2, 4 (bottom curve). The maximum is 
associated to the screening length Ls at the corresponding r]. 




Figure 3: L^TVcoshr] vs. rj. 



In the bulk of this paper, however, energetics are more interesting to us than 
screening lengths, so we turn to the computation of energy. The total energy is the 
sum of the various pieces, with the caveat that the energy of each string in the baryon 
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needs to be regulated by subtracting the energy of a free string stretching all the way 
from the boundary to the horizon. Before the regularization, the total energy is given 
by the (formal) expression 




Factors of two have been put in to take account of the two strings in each pair. The 
last piece is the energy of the D5-brane. The energy of the regulator quark has been 
calculated in the literature, we will use the expression (A. 10) from [4J. The idea is 
to cut off the energy integral written above at A instead of oo, subtract the quark 
energy integrated from the horizon to the cutoff, and then take the limit A to infinity 
after the subtraction. The result is. 



1 /"°° , / /l/^-p^cosh% \ 1 , , 



I r , ( ?/^-p^COsh% \ 1 J I - p4 cosh^ 77 

+ - / dy\ , - 1 + 1 - - + 



where we have used the variables 12.211 as well as Maldacena's relations and the geo- 
metrical definition of the Hawking temperature. 

With the values of L, a, (3 that we computed earlier, it is possible to make a 
numerical plot between L and E for various values of 7] and we present the results in 
the figure. This ties in as it should, with the similar plot in [5]. 

The computations of this section serve two purposes. Firstly and primarily for us, 
they give a context for the rest of the work on this paper. But in the process, they 
also provide a confirmation of the robustness of previous results on baryons. In [5], 
qualitatively identical results were observed, where instead of keeping the baryons on 
a circle, the angle at which the strings hit the D5-brane was fixed. This amounts to 
considering quark configurations that are squashed in the plasma wind. This makes 
the computations technically different, but still we see from a glance at the plots 
that the essential features are identical and therefore that the results are indeed very 
robust. In particular, the screening length formula L^Tv^cosh?] ~ const, holds in our 
case as well for large r], even though the value of the constant seems slightly different^]. 



^The slight variation in the constant is not bad - its value is known to depend on the details of 
the configuration, see e.g. figure 7, in [5j. 
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3 Dissociated Baryons 

Clearly there are many possible ways in which the baryon presented in the pre- 
vious section can dissociate. One possibility is to look at configurations where one 
quark has dissociated, while keeping as many of the remaining quarks as possible still 
on the circle. After the force balance conditions are imposed on the D5-brane, there 
are only a few interesting configurations which one can consider in this manner (we 
will look at them in Appendix C). But since the process of ripping a quark is funda- 
mentally dynamical, it is not clear to us that such configurations where the quarks 
are forced to be in a circle are the only ones preferred. Since it is difficult to come 
up with an unambiguous definition of what is "closest" and what is "farthest" to the 
original baryonic configuratioij^, we will attempt something more modest here. We 
will instead consider two baryonic configurations with identical linear dimensions but 
different orientations. We hope to make some comments about their relative ener- 
getics from this, but notice that we are steering clear of the question of what are the 
preferred dissociation products. 

The two configurations we consider are as follows. One is when the undissociated 

^Also, it is not clear how general our understanding would be even if we were to come up with 
such a notion for the four-quark baryon. 
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quarks are in a rigid identically spaced line parallel to the wind (i.e., along X3) (case I), 
and the other is when they are in the direction perpendicular to the wind (i,e., along 
Xi) (case II). We will also have a dissociated string in each case, which we will take to 
be well-separated from the rest of the baryon. Our aim is to look at two configurations 
that are identical, except for their orientation in the wind. Even though we will need 
to work much harder to get a fuller understanding of phenomenology, to get a hint 
about the orientation-dependence of the energy, this computation should be enough. 

We address both the cases separately. But before doing so, we notice a useful 
fact: the dissociated string that trails all the way to the horizon from the D5-brane 
is forced to lie along x^. This is demonstrated in Appendix B. 

3.1 Case I: Longitudinal Quarks 

We first consider the case where the undissociated quarks are along the direction of 
the wind. After the dissociation, the string corresponding to the dissociated quark is 
trailing, with one of its endpoints at the D5-brane vertex and the other at the horizon. 
The configuration we are considering is demonstrated in figure |5l and consists of three 
quarks in a straight-line. In case II, we will consider three quarks, again in a straight 
line, but perpendicular to the wind. The hope is that these two complimentary 
configurations will give us some idea about the energetics. Notice that because of the 
constraints arising from force balance at the D5-brane and because the trailing string 
is forced to lie along X3, many of the configurations are ruled out. An advantage of 
comparing the two cases we consider here is that they give a natural way to contrast 
between the cases: we can compare quark configurations of equal linear dimensions 
on either side. This is useful in getting some insight about how the energy of a given 
configuration (with fixed linear dimensions) changes with the orientation. 

We can compute the energy of this configuration using methods similar to those of 
last section. In what follows. Pi, ai are defined as in the last section (the indices in the 
present case denote the respective quark). In case I, the are identically zero. The 
origin of the coordinate system is taken to be at the D5 brane {xi{D5) = xs{D5) = 0). 
The numbering of the quarks are indicated in the figure. The trailing string solution 
forces the relation 

Pi = sinh?]. (3.1) 
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Figure 5: Case I: Remnant with quarks along the plasma wind. 



Using the results of Appendix A, we can write down the vertical force balance condi- 
tion at the D5-branc as 

Vl -p4 cosher/ / A r - — f - -Z\ 1 +p^cosh^?7 

and the horizontal balance as, 

/9i + /92 + /Sa = sinh r]). (3.3) 

In case I, we will look at configurations with the three undissociated quarks placed 
equidistantly on a line, along the wind. Thus we have the two equations below for 
the distance between quarks. 



^23 




These equations follow directly when integrating the equations of motion presented 
in Appendix B, for the specific case under consideration here. 
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An algorithm for solving the system is as follows. For each value of r] we do the 
following: 



• Pick a value of p. 

• For each value of p, scan /?i between -1 and 1. 

• For each value of Pi, solve (13.21) and (13. 3p simultaneously to obtain /32, /^s, paying 
attention to the signs using the geometry of the configuration. 

• Now we can compute L12 and L23 for each value of f3i. The value of f3i at 
which the two coincide is the correct value of jSi, and the corresponding L can 
be plotted as a function of p (for each given rj). Once we have fixed Pi, it is 
straightforward to fix P2, 03 using (13.21) and (13. 3p . 

This can be repeated for each value of 77. Once we have L and (3i as functions of p, 
we can compute the energy of the configuration, which is what we are really after. 
To compare the energies of cases I and II, we will not need to worry about the far- 
separated quark, so the regulated energy of the configuration shown in the figure 
takes the form: 



E = tTA 




y"^ — p"^ cosh 1] il,]^ l,vl~P'' cosh r] 



.(3.4) 



Here and in the next subsection, we have subtracted the energy of three regulator 
strings as opposed to four in the case of the undissociated baryon. We will present 
plot of L vs. p and E vs. L for both cases I and II together at the end of the next 
subsection. 



3.2 Case II: Transverse Quarks 

Case II corresponds to (undissociated) quarks aligned perpendicular to the wind as 
in figure [61 The Nambu-Goto Lagrangian for the string i that is stretched both in 
the X3 and Xi directions is given by 
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The energies, equations of motions etc. for quarks 1 and 3 in figure [H] are computed 
with this expression. 

Again, from the generic equations of motion written down in Appendix A, we find 
that the horizontal force balance conditions at the D5-brane enforces 

|a3| = |ai|=a, /^a = A = (3.6) 
/32 + 2/? = /34(=p^inhr/). (3.7) 

The vertical balance condition can be written as 

( 1 — cosh^ r]) l + p^cosh^T] 
(1 - P^) x/l -p4 cosher] 



2y'(l - p4cosh'r7)(l - p^ - (3^) - (1 - p^)a^ ^ {1 - p4cosh'r/)(l - p^ - /?|) 

The quarks 1, 2 and 3 all have the same X3 coordinate because of the configuration 
we have chosen. This means that 



p/?2 r 1 /?/4-p4 cosher] 



ttT a "(y^-p4)V y'-p'-A 

P^ r^^_ (y^-p^ cosher/) 



(?/4 - p4) ./ (y4 _ p4)(^4 _ p4 cogh2 ^) _ ^2(^4 _ ^4 ^Qsh^ ry) - a\y^ - p4) 



The spacing between the undissociated quarks (fixed to be L) is given by their xi 
coordinate, which can be obtained by integrating the equation of motion: 

L = ^ r dy ^ ( 3.10) 

"^1 ^J{y^ - p^){y^ - p^ cosh^ rj) - (3^{y^ - cosh^ r]) - a^{y^ - p^) 

These equations can again be solved numerically for any value of r]: 

• Pick a value of p. 

• For each value of p, scan [^2 between and 1. 

• For each value of /32, solve (13.71) and (13. 8p simultaneously to obtain a, 13. 
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• Now we can evaluate the RHS and LHS of (I3.9p for each value of /?2- The value 
of (32 at which the two coincide is the correct value of (32- Once we have fixed 
P2, it is straightforward to fix a, /5 using (13 .7^ and flS.Sp . The corresponding L 
can be plotted as a function of p (for each given rj) using (13. 101) . 

The regulated energy of the dissociated case II configuration can be calculated 
similar to the previous cases (here again, we subtract three quarks as in case I) and 
the result is 



E = tVX 



dy 

P Ji 



{y'^ — cosh^ rj) 



1 + 



{y^ - p4 cosh^ r]){y^ - p^ - p^) - a^{y^ - p^) 

y'^ — p^ cosh^ r] \ 1 a/1 — p"' cosh^ 77 



y 



-1+1 



P 



+ 



2p 



.(3.11) 



3.3 Plots 



Using the expressions from the previous subsections and the numerical simulations 
resulting from them, we can make comparisons between the energies of the two disso- 
ciated configurations. The curves are similar in both cases, except that the longitudi- 
nal quarks are at a higher energy when the windspeed is non-zero. The energy plots 
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0.4 - 




Figure 7: The vs. L plots of the transverse and longitudinal configurations at 
1] = 0. Even though the numerics that gives rise to the curves is different, the two 
are exactly on top of each other. 




0.0045 



Figure 8: The vs. L plots for the two cases at = 6. The longitudinal curve 
is the one at higher energy. The qualitative features are identical at other values of 
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0.005 r 




Figure 9: LT vs. p for r] = Q for the two cases. The transverse case (II) is marginally 
lower, as expected from the intuition of figure 4 in [5]. 



0.4 r 




0.00 0.02 004 006 0.08 0.10 n.ono 



(a) r;=0. (b) 77 = 6 

Figure 10: Plots of vs. LT for baryons (lower curves) and dissociated configura- 
tions (upper curves). For less clutter we only show the transverse case. 

should be identical when the windspeed goes to zero since case I and II are identical 
in this case, so this can be used as a consistency check of our numerics. This is indeed 
what we find as clear from the r/ = plot. We show the plots of E vs. L for two 
representative cases (r/ = and = 6), in their interesting regimes of parameters. 
The fact that longitudinal configurations are at a higher energy than transverse ones 
has previously been observed in the case of mesonic configurations in [1] (see right 
panel of figure 6 in |1] ) . 
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Another observation one can make from the plots is that the energies of the 
dissociated configurations are always above those of the corresponding undissociated 
configurations for sufficiently large values of rj, and this is shown in Figure llObI This 
result is reasonable. The zero of the energy at each rj is set by the (colored) trailing 
string at that rj. We also know from a previous section (see also [H]) that baryons have 
negative energy with respect to this datum, for sufficiently large windspeeds. So the 
configurations we are investigating here, which are morally between baryons and free 
strings, should naturally have intermediate energies. For small enough windspeed, 
there is some changes in these comparative plots, and we present them in Figure llOai 
The cusp region is above zero energy even for the undissociated baryon in this case, 
as expected (see also [5]). 

The structure of the energy curves for the dissociated cases is roughly 
similar to that of the undissociated case, but in 
the plots here, we have zoomed in on the fine 
structure. The one significant difference with the 
baryonic case is that instead of a cusp, for non- 
zero boost, here we see a loop. We have done these 
numerical simulations for various configurations 
and various resolutions (see also Appendix C), so 
it seems very unlikely that this loop is a result 
of some systematic numerical error. The loops 
robust feature of all dragging configu- 
rations. In particular, they vanish for color non- 
singlet configurations only when the windspeed is 
zero and there is no drag (see Figure [T]). The tech- 
nical reason why the loops arise in E vs. L plots 
is not hard to see. The schematic plots of E vs. p 
and L vs. p for generic configurations and generic 
rj are given in the figure to the right. They both 
have one maximum, and for the case of non-drag 
configurations, the peaks in both E and L happen 
for the same value of p. But for dragging config- 
urations, the peaks in E are displaced to the left 
with respect to the peak in L, and therein lie the 
origins of the loop. It would be nice to understand the physical origins of this shift. 




Figure 11: Typical shapes of E 
and L vs. p for any rj and for 
any configuration. The axes have 
been re-scaled to fit the curves on 
the same plot, and the specific 
values are not to be paid atten- 
tion to. 
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better. 

We suspect that these loops can be used as a measure of the drag of the configu- 
ration, because their areas seem to vary depending on the configurations l|. The fact 
that the loops have non-zero area for dragging configurations, and the fact that the 
area of the loop in the vs. LT plot has dimensions of energy in string units, 
suggests that this could perhaps be used as a measure of the energy loss. It would 
certainly be very interesting to study this further, and we plan to come back to this 
in the future. 



4 Discussions and Loose Ends 

We have already reported the main conclusions in the introduction, so here we 
will merely make some comments. 

There are some natural extensions to the work done here. We have considered the 
special case of A^^c = 4 quarks, and found that we can reproduce and extend the results 
in the literature that deal with color singlet configurations. Considering a uniform 
distribution of quarks along the circle where the large limit is more systematic 
would be a natural next step. This will be more complicated because the strings 
are now not along the coordinate axes. More generic gauge theories, configurations 
other than circular baryons, more generic wind directions etc. are all possible lines 
of investigation, especially in understanding the possibility of extracting universal 
predictions valid across conformal/confining gauge theories. 

In this paper, we have compared dissociated configurations which capture some 
aspects of the energetics of baryon dissociation. But it should be emphasized that the 
results of this paper are tentative. It would be nice to investigate other configurations, 
especially those with more quarks. Another interesting possibility is to do a scan of 
the various dissociation configurations at various angles and lengths. In particular, 
investigating other configurations which are identical except for their orientations 
would be useful. 

Our original configuration of four quarks was specifically chosen to simplify the 
calculations. This is a perfectly good starting point, but it also suffers from the draw- 
back that we cannot investigate many interesting configurations. A good example is 



^For the cases considered in this section, this effect is not pronounced. 
^Not all such configurations are kinematically allowed. 
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the transverse and longitudinal configurations (analogous to the two configurations 
we investigated in section 3) after two quarks have dissociated. Both these configu- 
rations are not allowed because the force balance conditions on the D5-brane are too 
restrictive. 

To get a full and unambiguous understanding of the possible dissociation patterns, 
we need a more exhaustive study of the various hints we have found in this paper. 
We have tried to stay close to the configurations which are most easily tractable. Our 
aim here has been to set up some of the framework and leave the more thorough work 
to a more elaborate future project, perhaps with more man-power. Some of these 
questions are currently under investigation. 
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Appendix 

A. Equations of Motion 

The equations governing the various static baryon-like configurations is obtained by 
varying the total action with respect to the string coordinates and with respect to 
the position of the D5-brane. Since the boundaries of the strings (at the D5-brane) 
are also supposed to be varying, one ends up getting a bit more than the usual Euler- 
Lagrange equations. The basic ideas are presented in [5], but we have chosen to 
redo it here for the case when there are dragging strings to clarify the origin of some 
negative signs which turn out to be crucial in this work. 



22 



The action for the system takes the general form 



E 

ie{up} 



ne{down} 



Cn{Xi„,r)dr + SD5{re), (A.l) 



where we have denoted all the functional dependencies on the relevant variables. Xi^rn 
stands for the the i-th coordinate of string m. Primes are, again, with respect to r. 
The summation over the strings will be written explicitly in what follows, but the 
summation over i, should be understood from the context. 

First wc vary with respect to the x's and get the individual equations of motion 
for the various strings. But since the boundaries arc also allowed to vary, we also get 
boundary equations of motion which give further constraints on the configurations. 
Setting SS = (for variations of Xi), we end up with 

dC 







E 

me{up} 



dr 



= E 

me{up} 

d ( dCr, 



°^ dC 



ne{down} 



dr 



dx 



-5x',. 



(A.2) 



i,n 



dr V 9a;!!- , 



-Sxi 



d ( dCr 



dr V dx', , 



5xi . 



+ E 

nS{down} 



dr 



m 



n 



We have written 5x'^ ,^ as d{5xi^rn) / dr and done an integration by parts as usual. Since 
the variations Sxi are arbitrary in the bulk (of the string), the second term in each 
piece has to be zero and gives rise to the standard Euler-Lagrange equations. In the 
present case they take the form: 



const. = Kj . 



What remains in the variation is 

dCrr. 



dx',. 



const. = Ki , 



(A.3) 
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+ E 

ne{down} 
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dr V dx'i , 



■5xi 



= - E 



mG{up} 



dx',, 



-Sxi 



+ 



E 



dC 



dx'- 

ne{dowii} 



— SXi- 



(A.4) 



The strings are fixed at r = rg and r = cx), so the terms from those ends of the integral 
don't contribute. Also, at r = r^, the strings are allowed to move, but only under 
the constraint that they are all still attached at the baryon vertex. This means that 



6xi, 



m\re 



5xi^n\re = Sxi\r^. Another input comes from the Euler-Lagrange equations 
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above. They say, in particular, that = Ki. Putting this all together, we end up 

with 

Ki^m- Yl ^i^n = 0, (A.5) 
m£{up} nG{down} 

which gives rise to three equations, one for each i. 

Now we turn to variations in rg. The idea here is that we vary the boundary re, 
an look for variations in x which result from considering the extrema of the action 
with this new boundary. In other words, in Figure [T21 both x and x are solutions of 
the Euler-Lagrange equations with their respective boundaries. 




Figure 12: The classical solutions corresponding to different boundaries. 



One way to handle this shift in boundary is to think of r as a map from u G 
[0, 1] to the integration range (i.e., write r as r(-u)). Then we can think of u as the 
time variable, with fixed boundaries and 1, and r will be just another coordinate. 
The changes in the boundary can be re-interpreted now as variations in r, at the 
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fixed boundaries of tlie w-interval. Tfius we liave translated a moving boundary to a 
boundary term: sometliing easily liandled using standard variational approacfies (See 
e.g., cfiapter 1 of p31 for a clean discussion of variational methods with boundary 
terms.). This met ho co works, but unfortunately our action involves different pieces 
which contain different ranges for r: the upward strings go from r^, to oo, while the 
downward strings are from tq to Tg. This makes this approach somewhat complicated, 
so we will follow another path which is much more direct and intuitive. 

Let us consider what happens to the classical solutions of a system with the action 

So = r L{x{r),x\r))dr, (A.6) 



under variations of the boundary. We are interested in finding new equations of 
"motion" by setting 

pre+Sre pre 

Q = S^,^So= L {x, x') dr- L {x, x') dr (A.7) 



The important thing to note here is that x is supposed to be the classical solution 
corresponding to variations with the boundary fixed at + Sr^, just the 
classical solution with boundary fixed at rg. We have denoted x{r) — x{r) = 6x{r). 
Upto first order in small quantities, we can write, 

rre 

0= {L{x,x') - L{x,x'))dr + L{x,x')\rjre. (A.8) 
Jo 

By the usual tricks, the first piece can be massaged into the form 



re / QJ^ QJ^ 

dr I — 5x + — — 5x \ = I dr 
ox ox' 



\ dx dr dx' J dr \ dx' 



The Euler-Legrange term vanishes because x is a classical path, and the other piece 
gets integrated and receives contributions only from the boundary. Since we imagine 
that the boundary at zero is held fixed, we end up with 



dx' 



+ L{x,x%Jre. (A.9) 



An important input at this stage comes from the fact that 6x{r) is the variation at 
fixed r between classical solutions with different boundaries. Notice that by definition 



^Incidentally, this approach seems like a pretty general way to do variational mechanics. 
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STx{rf.) = Sx\r^ + §7^^"!^ is zero: at the new boundary r^, + dr^, the new classical 
solution X is supposed to have vanishing variations as well. So we get the final result, 







dx' 



(A.IO) 



This was done for the case when the moving boundary was at the upper end of 
integration. For quarks hanging from the boundary of AdS, there will be an overall 
sign. Putting all these ingredients together, we finally get 



mG{up} ^ '■'^ ^ ^-^^ nG{down} ^ 



+ ^ = 0(A.ll) 



as the Te-equation of motion for the various quarks and the D5-brane. 



B. Trailing String 

In order to get a handle on the various dissociated baryonic configurations, we need 
to understand the trailing string solution that extends from the D5-brane to the black 
hole horizon. We will show in this appendix that the most general string of this form 
in the Xi — x^ — r space lies along the wind direction {x^ — r plane). This intuitively 
natural conclusion is important because (due to the force balance conditions at the 
D5 brane) it considerably restricts the resultant dissociation configurations that are 
allowed. 

The action for the most general static string stretched between the D5 and the 
horizon is. 



s= ^ 



The equations of motion for the two components are 



By defining 



= — , (B.13) 

x' = ^^-^ (B 14) 

^ fr y/r^fA - R^KIA - R^Kff 



r a j3 

— , A* = ^, = ^ and a;i,3 = — 2:1,3, (B.15) 
ro p^ ro 
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we can rewrite the above equations as 



dzi \ ^ 2 1 



, , ^ , (B.16) 

fdz,y ^ (t^-cosh^)^ \ 

We would like to have simultaneous solutions of these two equations. A real solution 
can clearly only exist if 

(t^-cosh2r7)(t^-l-z/2)-/i2(t4_i) >0. (B.18) 

(We will be sloppy about distinguishing > and > in what follows, we will consider the 
boundaries explicitly.) The variable t ranges from t = 1 at r = tq (horizon) to t = 1/p 
at r = Te (D5 brane). We are interested in the case where '\-/p> V cosh r] (See, for 
example, the screening length plots.). So for the entire range 1 < t < \/ cosh rj + e, 
with e small enough0 and positive, the above inequality should hold for the solution 
to make sense. Doing a variable redefinition, this means that 

+ a;(sinh^ rj — jjL^ — z/^) — jj? sinh^ ^ > 0, (B.19) 

should hold for — sinh^ rj < x < 5 with small, positive 5. The above inequality is 
satisfied as long as x is either less or more than both the roots of the quadratic!"] 
(the discriminant is positive, so the roots are never complex). If one of these ranges 
overlaps with — sinh^ r] < x < 5, then we have a solution. This happens iff 



— (sinh rj — ji^ — z/^) + ^/(sinh rj — jj? — u'^Y + Ajj? sinh rj 
- sinh^ rj > ^- , (B.20) 



or 



— (sinh rj — — u ) — \ (sinh rj — — u^y + 4/1^ sinh rj 

6< ^ . (B.21) 

2 ^ ' 

The former inequality implies 

z/2 sinh^ 7] < 0. (B.22) 



^''The integration limits are between 1 and i. We know that ^ is bigger than ^/coshr], even though 
we don't know how much bigger. So at least for sufficiently small positive e, we can claim that the 
inequality should hold. 

^^As opposed to being lesser than the bigger root and bigger than the smaller root. 
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This is clearly impossible, but there remains the possibility of the limiting case = 0, 
when the inequality is saturated. But this case is clearly pathological because it is 
easy to see that the original inequality (IB.ISI) is violated for values of t that lie close 



enough to 1. The remaining possibility is the second inequality (]B.21[1 . This can be 
rewritten as 

+ ^(sinh^ r/ - /i^ - z/2) - sinh^ r] > 0, (B.23) 

and it needs to be satisfied if we tune 6 small enough. This means that the only 
possibility is the limiting case /i = 0, where the inequality is saturated. But this is 
precisely the case when the trailing string is purely along xs, which gives rise to the 
solution considered in appendix A of [1]. 



C. Other Configurations 

In this appendix we will look at some dissociated configurations to illustrate the 
fact that the energy plots can significantly differ from the ones considered before, 
depending on the configuration of quarks. The configurations we consider here are 
the minimal configurations allowed, if one stipulates that after a quark has dissociated, 
the remaining quarks in the baryon still remain on the circle. Of course, since we 
do not understand the dynamics, this is an ad hoc assumption. The purpose of this 
section is to give a flavor of the various dissociations patterns allowed. 

0.006 r- 
0.005 
0.004 
0.003 
0.002 
0.001 
0.000 

0.00 0.01 0.02 0.03 0.04 0.05 0.06 

Figure 13: LT vs. p for the various cases. The top curve is for the undissociated 
baryon. The color code for the rest is as follows. Violet: case B, blue: case A, green: 
case I (from section 3), brown: case II. All the plots in this appendix are for rj = 6. 




28 




0.0035 0.0040 0.0045 



Figure 14: vs. L at r/ = 6. The color code is the same as before. 




The configurations we consider are easily described by the condition that the 
quarks remain on the circle. The force balance conditions for our simple system are 
stringent enough that if the dissociated quark is one of the transverse ones, then not 
all of the remaining three quarks can remain in a circle. So we will look at the two 
cases where one of the quarks in the longitudinal direction is the dissociated on^. 

There are a few configurations one can consider in order to get an idea about transverse quark 
dissociations. One of them was already considered in section 3. Many other possibihties are ruled 
out or become uninteresting because of the constraints on the configuration. For example, the case 
where all three undissociated quarks are at the same location at the boundary, is fully fixed by force 
balance and symmetry. Also, it does not have the analog of a screening length, so it is not very 
interesting for our purposes. 
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In the (xijXs) plane, the two cases we consider will be defined by the quark positions 

Case A: (1, 0), (-1, 0), (0, 1), (C.24) 
CaseB: (1, 0), (-1, 0), (0, -1). (C.25) 

For example, in this notation, the longitudinal case (case I) considered in the main 
text will be defined by the quark locations (0, 1), (0, 0), (0, —1). 

We will not present the details of the computations because the general formahsm 
is the same as in the examples we considered in the main body of the paper. The one 
slight subtlety is that to make sure that the boundary quarks are arranged a circle, 
the relations we need to impose are, 

4? - 4i °^ - 4'. (C.26) 

for case A, and 

for case B. The L^^ are measured from the D5-brane, and the superscripts denote 
the relevant quark (notice the 1^3 symmetry). We present here the plots of the 
energies, screening lengths, together with those of the cases considered in section 3. 
The comparisons with the undissociated baryons are also presented. 
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